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ABSTRACT 

This paper presents the complete forward and inverse kinematics solutions for 
control of the three degree-of-freedom (DOF) double octahedral variable geometry truss 
(VGT) module as an extensible gimbal. A VGT is a trass structure partially comprised of 
linearly actuated members. A VGT can be used as joints in a large, lightweight, high load- 
bearing manipulator for earth- and space-based remote operations, plus industrial 
applications . The results have been used to control the NASA VGT hardware as an 
extensible gimbal, demonstrating the capability of this device to be a joint in a VGT-based 
manipulator. This work is an integral part of a VGT-based manipulator design, 
simulation, and control tool. 


1. INTRODUCTION 


Variable geometry trusses (VGTs) are a dass of in-paralld actuated devices. In- 
parallel actuated manipulators have several load paths from the end-effector to the fixed 
link, as opposed to serial manipulators which have only one. They are truss structures 
with linear actuators for some members. VGTs retain the desired stiffness and load- 
bearing characteristics of trusses, while providing the capability to articulate. VGT 
manipulators are lightweight, have long reach, and are high load bearing compared to 
equivalent serial manipulators. All actuators are linear, which are preferable to rotary 
joints. VGT manipulators can be constructed in a modular fashion, allowing many 
redundant degrees of freedom (DOF). This redundancy can be used to optimize 
performance, such as providing snake-like motion to avoid obstacles, avoiding 
singularities, and avoiding joint limits. The VGT has an open structure allowing routing of 
cables and hoses. When properly designed, all VGT members are loaded axially, thus 
increasing stiffness and load bearing capability. 

The double octahedral VGT (Fig. 1) was developed at NASA Langley in the 80’s 
(Rhodes and Mikulas, 1985). It was intended to be used in large deployable space 
structures. VGTs were then pursued as active structural elements to damp vibrations of 
large space structures (Wynn and Robertshaw, 1991). Recently, researchers have 
considered using VGTs as robotic manipulators (Mill sap, et. al., 1992). These authors 
propose a large-scale VGT-based manipulator to clean up large waste storage tanks. The 
Department of Energy (DOE) determined that the state of the art in large manipulators 
was not sufficient to efficiently handle this problem (Mill sap, et. al., 1992). Another DOE 
application for VGT technology is as a trunk for carrying two dexterous serial 
manipulators on a decontamination and decommissioning vehicle for cleaning up 
contaminated nuclear facilities, under development at Oak Ridge National Laboratories. 

The VGT shown in Fig. 1 has three DOF. It can be used as a joint in a modular 
manipulator constructed of articulating joints and static truss sections. In this paradigm, it 
is convenient to control the VGT module as a gimbal, with a, (3 rotations about 
perpendicular axes. In addition, accordion-like extension of the VGT joint is possible. 
This extension can be fixed at a nominal middle value (i.e. constrain the three DOF to 
two). However, to obtain the most capability, all three DOF should be coordinated. The 
purpose of the current paper is to present the complete forward and inverse position 
kinematic solutions to control the double octahedral VGT (hereafter referred to as VGT) 
as an extensible gimbal. 


2 



The groundbreaking theoretical kinematics work for this VGT was performed at 
Virginia Polytechnic Institute & State University. Padmanabhan (1992) presents a general 
study of VGTs (called isostatic frameworks) in his PhD dissertation. The double 
octahedral VGT is a case study in this work. Padmanabhan, Tidwell, Salerno, and 
Reinholtz (1992) present the concept and general solutions for using the VGT as a gimbal. 
Padmanabhan, Aran, and Reinholtz (1990) present a general closed-form solution to the 
VGT inverse position problem. 

The primary contribution of the present paper is that it summarizes all kinematics 
equations necessary to model the VGT as an extensible gimbal, including full details. 
Original work reported includes formalization of the subproblems for forward and inverse 
position kinematics, Stewart platform analogy for the forward solution, closed-form 
expressions for the improvement vector at each Newton-Raphson iteration for the forward 
solution, VGT mobility calculation, VGT node locations calculation, and presentation of 
the complete forward transformation for the extensible gimbal. 

The current paper is organized as follows. Section 2 summarizes the symbols 
used throughout the paper. Section 3 presents the paradigm for controlling the VGT as an 
extensible gimbal. Section 4 presents the complete forward and inverse position 
kinematics solutions for the VGT, including the extensible gimbal equations. Section 5 
gives examples to demonstrate the various solutions. Appendix A presents mobility 
calculations for the VGT, Appendix B presents fixed VGT geometry, Appendix C applies 
the Newton-Raphson method to the VGT forward kinematics solution, Appendix D 
provides extensible gimbal details, and Appendix E presents a method to calculate the 
VGT node locations for finite element modeling and general simulation. 
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NASA VGT Hardware 
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2. SYMBOLS 


a,;/ = 1,2,3 
A, B, C, D, E 

B IR,B 2 r>&sr 

{B}, {P} 

Cj, Sj 

A » B>i > * * = 1» 2, 3 
Ei,F„G,;i = 1,2,3 
L 
Lo 

L 1 ,^ 2, ^3 
Emin ' E MAX 
n = X P 

{n x n y it;} 7 

n x 

N,N it N 2 ,N 3 

B ap<h,<h t <h 

Oi y ,Oj z ;i = 1,2,3 


P 

Po 

b 

I 


Qi-Qi'Qz’Qi 


QQj 

QlR’QlR’QhR 

r 

s 

[fr] 

f>i=(ui, <4, 

x* _ 

%B’Yb’Zb 

Xp,Yp,Zp 

a,P 

(e.), 


'U 


lull 



Vectors locating = 1,2,3 in actuation plane 

Simplified forward kinematics coefficients 
Base nodes 

Top platform nodes 

Base and platform coordinate frames 
cos0 ; , sin0, 

Forward kinematics coefficients 
Inverse kinematics coefficients 
Cross longeron length 
Fixed batten length 
VGT input leg lengths 
Minimum, maximum actuator length 
Top platform normal 
Components of n 

1 + n x 

Cross longeron face normal length 

Stewart platform revolute locations 
Components of 0 lr 0 2 ,0 3 

Length from {BJ to {P} 

Difference of p and joint offset S 

Cross longeron face unit directions 
Lower actuation plane nodes 
Length from Q to Qj 
Upper actuation plane nodes 
Gimbal extension parameter 
VGT mid-plane joint offset 

Homogeneous transformation matrix, {P} in {B} 

Mid-plane normal 
Base normal 
Unit vectors of {B} 

Unit vectors of {P} 

Gimbal angles 

Angles of projections of q lf q 2 ,q 3 onto base 
Cross longeron face angles 

Two inverse position solutions to 0, 

Euclidean norm of vector x 
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3. DOUBLE OCTAHEDRAL VGT AS AN EXTENSIBLE GIMBAL 


The NASA VGT module shown in Fig. 1 can be controlled as a two-axis gimbal 
(universal joint) with extension capability. Figure 2 shows the kinematic model for the 
VGT. As shown in Appendix A, the VGT has three degrees of freedom (DOF): there are 
three independently controlled linear actuators To provide three DOF, the 

design of joints connecting links at the VGT mid-plane is critical. Rhodes and Mikulas 
(1985) present the first VGT joint designs based on slotted spherical joints. Tidwell 
(1989) presents another design based on universal joints. Padmanabhan, et. al., (1992) 
present a good overview of VGT joint design. The NASA VGT hardware has mid-plane 
joints constructed with universal joints, but the details differ from Tidwell (1989). A 
promising alternative spherical joint design has recendy been proposed by Hamlin and 
Sanderson (1994). 

Since the VGT has three DOF, three output parameters can be controlled. There 
are five possible Cartesian output parameters to control ( x, y, and z translations of a 
point on the platform with respect to the base, and pitch and yaw rotations about the Z P 
and Y P axes, respectively). The roll rotation about X P is not possible with the VGT. For 
motion control, one possibility is to control the translation of a point on {P} with respect 
to the origin of { B }, as presented in Padmanabhan, et. al., (1990). With this scenario, the 
rotations cannot be controlled. 

In this paper, the three output parameters are the pitch and yaw angles, plus the 
length of the joint. The VGT in Fig. 2 is modeled as the two-axis gimbal shown in Fig. 3. 
The variable length of the virtual gimbal joint, r, is the third freedom. This extension is 
not available in conventional gimbal designs which provide only the pitch and yaw angles. 
In addition, the stiffness and load-bearing characteristics of the VGT are superior to 
conventional gimbals due to the parallel truss nature of VGTs compared to the serial 
gearbox design of conventional gimbals. The intersection point of the two gimbal axes is 
fixed with respect to the VGT model. The lower and upper gimbal joint lengths (r in Fig. 
3) are equal due to symmetry. 

The three degree-of-freedom VGT can be used as a joint in a VGT-based 
manipulator. Such a manipulator could be controlled as a series of pitch-yaw joints with 
extension. The extension r can be constrained to be a constant, intermediate value for 
simpler manipulator kinematics. This case effectively maps the three degrees of freedom 
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to two output parameters (pitch and yaw). However, better performance may be possible 
if the extension r is also actively controlled. 



Figure 2 

Double Octahedral VGT Kinematic Diagram 




Figure 3 

Extensible Gimbal Kinematic Diagram 
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4. DOUBLE OCTAHEDRAL VGT POSITION KINEMATICS 


This section formulates and solves the complete forward and inverse position 
equations for the VGT. The VGT input parameters are Lj.L2.L3 and the output 
parameters are a,p,r. The forward problem calculates a ,p,r given Lj.L i,L$, while the 
inverse problem calculates Lj.L2.L3 given a,(5,r, as illustrated in Fig. 4. The forward 
solution is useful for simulation and the inverse solution is used primarily for control. 
Both are used for VGT design and analysis. 



Figure 4 

Position Kinematics of the VGT as an Extensible Gimbal 

Figure 5 shows intermediate details of Fig. 4. Each block in Fig. 5 is described in 
the following subsections. The VGT Forward A solution of Section 4.1.1 calculates the 
cross longeron face angles 0j.02.03 given the input parameters Lj.L2.L3. The VGT 

Forward B solution in Section 4.1.2 calculates the pointing vector h and the gimbal 
extension r given 0 j, 02,03. To finish the forward solution, the Gimbal Inverse solution in 
Section 4.3.2 calculates the gimbal angles a.P given h ; the gimbal extension r is 
unchanged. This solution is inverse because a.P are gimbal inputs while h is output. The 
VGT inverse solutions are accomplished using the same intermediate variables, in reverse 
order. The Gimbal Forward solution is presented in Section 4.3.1, VGT Inverse A in 
Section 4.2.1, and VGT Inverse B in Section 4.2.2. 

Both the VGT forward and inverse solutions exploit the symmetry inherent in the 
design of Fig. 2. As explained in Padmanabhan, et. al., (1992), the actuator plane 
Lj — L2 - L3 is a plane of symmetry. Due to the specific VGT geometry, the motion of the 
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top platform is a mirror image of the base, reflected about this plane of symmetry. 
Likewise each VGT longeron has a symmetric partner about this plane. 



Lj L 2 L 3 

VGT 

Inverse 
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I Extensible 
Gimbal 
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Figure 5 

Position Kinematics of the VGT as an Extensible Gimbal: 
Intermediate Variables 


There are five length parameters for this VGT design (see Fig. 2): (1) Six bottom 
and top fixed battens Lq. The geometry is shown in Figs. B.l and B.2; (2) Twelve fixed 

cross longerons L. The geometry for the cross longeron faces is shown in Figs. B.3 and 
B.4; (3) Minimum actuator extension L MIN ; (4) Maximum actuator extension L MAX . 

The three actuators are identical; (5) Joint offset S separating pairs of cross longeron 
faces about the plane of symmetry. If a different parameter set was chosen for a VGT 
design, die desirable symmetry condition might not result. 


4.1 VGT Forward Position Solution 

This section presents the two VGT forward position solutions; given 

calculate h,r. The gimbal inverse solution (Section 4.3.2) is then required to complete the 
VGT forward solution. As shown in Fig. 5, the VGT forward solution is divided into 
solutions A and B, separated by intermediate variables 0j , 0 2 , 63 . 

4.1.1 VGT Forward A. In this section, the solution for cross longeron face 
angles 0 1 ,0 2 ,03 is obtained, given input parameters 3. Only the half of the VGT 

under the plane of symmetry is solved. Then the position of the top half is determined by 
symmetry. This solution was presented by Padmanabhan, et. al., (1992), where three 
RSSR mechanisms (R-revolute joint; S-spherical joint) 0iQiQ 2 0 2 , G 2 £? 2 £? 3 03 . 0^Q?,Q\0\ 
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were identified and solved simultaneously (see Fig. 6b). An equivalent solution recognizes 
that the VGT configuration under the plane of symmetry is a Stewart's platform, with 
different actuation: the cross longerons are fixed and the platform lengths are actuated. 
The forward position solution for the Stewart platform is presented in Williams (1992). 
This solution is adapted for the VGT below. 


T S/2 

Q 3 Qi 


Figure 6a 

Stewart Platform Analogy : Offset Detail 

Qi 



Bi 

Figure 6b 

Stewart Platform Analogy 

Figure 6 b shows the Stewart platform analogy for the lower VGT half. The VGT 
joint offset S does not affect this solution because the plane connecting points £?i-£? 2*£?3 is 
parallel to the plane of symmetry. Therefore, virtual links Lj,L connect points 
0i, 0 2 , Q 3 to form the Stewart platform; see Figs. 6 a and 6 b. Figure 6 b shows the 
definition of the cross longeron face angles 0 i, 0 2 , 03 , measured from the base plane to the 
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virtual links OQ , about the base fixed batten directions. In Fig. 6b, the solid links are 

virtual, while the dotted links are actual VGT links. All links are fixed except for 
which are given. 

Three scalar constraint equations are written, stating that L y , 1^,1^ are the lengths 
separating points Qi,Qi,Qy 


ii=Wi=\ B Qi- B Qii 

i|=a&=| B a- fl af u> 

l l=aG=|®a- e e 3 f 

where: 

B 0i= B Oi +Nj B qi . (2) 

The fixed locations B O t and lengths N t are given in Appendix B. The variable unit vector 
directions B qj pointing along OQ are functions of 0 1( 02,03, given in Eq. 3 (s t = sin0,- and 
c, = COS0,). 


\ = i 


■» i 

-c f c«t>, 


(3) 


From Fig. B.l, <►, = 180°, $ 2 = 300°,<t> 3 = 60°. Therefore, 
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< 7 1 = 


s 2 


*3 

1 

B 

1 

1 

1 

<73 =< 

c 3 

2 3 

V3 


V3 

IT C2 J 


c 3 

2 3 J 


II 

<5? 

Nsi 

Ojy + NCl 

, b Qi = < 

Ns 2 

°2y~\ c 2 

v 

to 

Kji 

11 

A. 

Ns 3 
- N 

Oly-~ C 3 ' 


0 


o + Sn 
r 2j+ 2 C2 i 


n V3 N 
[° 3 *“ 2 3 1 


(4) 


( 5 ) 
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Upon substitution, the three constraint equations become: 


D l c \ + D 2 c 2 + D 3 c l c 2 + A^l$2 +D 5 =0 

E x c 2 +E 2 c 3 +E 3 c 2 c 3 + E 4 s 2 s 3 +E s = 0 
F l c 3 + ^2 C 1 + ^3 C 3 C 1 + E 4^3^l +F 5 =0 


where: 

D l = 2N(0 ly -0 2y ) 

D 2 = - 0 2y + V30 2z ) 

d 3 = n 2 

D 4 =~2N 2 

D 5 =(o ly -0 2y f +01 + 2N 2 -! 2 

Ei = -0 3y )+ V3 (0 2z - 0 3z )) 

E 2 = N {{o 2y - o 3y ) + V3(0 2z -0 3z )) 

e 3 = n 2 

E 4 = ~2N 2 

E 5 = ( 0 2y - 0 3y f +{0 2z - 0 3z f + 2N 2 -I 2 
F i = N {0 iy -0 3y -y/30 3z ) 

F 2 = 2N(0 ly -0 3y ) 

F 3 = N 2 
F 4 = -2 N 2 

F 5 =(o ly -0 3y f +0 3z +2 N 2 -I 2 


The three constraint equations simplify to: 


Ci + c 2 + Ac x c 2 - 2As x s 2 + B = 0 
c 2 + c 3 + Ac 2 c 3 - 2As 2 s 3 +C = 0 
c 3 + Ci + j4c 3 Ci ~ 2As 3 s x + D = 0 



where: 


A iSn 
~ A " 3A 

*=t=£(^) 

c=|=^( E - 4 ) 

£> = ^- = -^■(£-4) 

A w 2 v ' 

E = -J-+2N 2 
4 


Equations 7 are three coupled non-linear equations in transcendentals of the 
unknowns 0i, 02,03- Due to the symmetric VGT geometry, the coefficients of Eqs. 7 are 

simple compared to Eqs. 6. However, the basic structure of the solution is unchanged: 
there are 16 theoretical solutions (Nanua, et. al., 1990). Actually, there are 8 potentially 
useful solutions; the remaining 8 solutions are symmetric about the base plane. 

For most simulation and control purposes, only one solution (closest to the 
previous configuration) is required. Therefore, in this paper, Eqs. 7 are solved numerically 
using the Newton-Raphson method. Appendix C presents adaptation of the Newton- 
Raphson method to solve for 0j ,02,03- 


4.1.2 VGT Forward B. The pointing vector h and the gimbal extension r are 
calculated in this section given the cross longeron face angles 0j,02,03, where ft is the 

normal to the VGT top platform, at its centroid. Symmetry is exploited to obtain the 
solution. This general solution was presented by Padmanabhan, et. al., (1992). Figure 7 
gives the vector diagram for this solution. 

Given 0!,0 2 ,03, B Qj are calculated using Eq. 5. Then the unit vector directed 
from the origin of {B} to the origin of {P} can be found with Eq. 8. This unit vector U l is 
normal to the actuation plane of symmetry. The total length from the base to the top 
platform along this vector is p. 


A = 


(*g 2 - B giM*g3-*gi) 

b Q 2 - b Qi)4 b Q3- b Qi 


I A, 

Ay 
I Az . 


(8) 
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The projection of r along U l is this yields the following condition: 


r = — P_ 


2X B U l 


(9) 


where X g = {1 0 0} r , thus: 


r = 


2f/i 


lx 


( 10 ) 


To calculate r, we need p. The projection of B Q t onto f/j is: 


A) fl 


Qi-u !• 


( 11 ) 
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Since is normal to the mid-plane, the joint offset S is parallel to it. Therefore, the 
distance to the mid-plane along U l is: 


Z = £o+^ 

2 2 2 ' 


( 12 ) 


Substituting Eq. 1 1 into Eq. 12 and simplifying yields Eq. 13. Using this expression for p, 
r is determined from Eq. 10. 


p=2 B Q v U l +S (13) 

To determine n , the following vector loop closure equation is written. 

rX B +rn = pU l (14) 

Solving, 

« = ^U 1 -X B . (15) 

r 


The expression for it in Eq. 15 is a unit vector, as proven below (making use of Eq. 10). 


h = pWu 


U\—X B = 2U\ r U) -X R = 






2Ut-\ 
2U lx U ly 
2 U lx U lz 


(16) 


W = ufful* K + vh - 1) + 1 = ■ + 1 - 1 (H) 

For each of the sixteen possible 0 1 ,02,03 solutions from Section 4.1.1, there is one h,r 
solution. 


4.2 V(»T Inverse Position Solution 

This section presents the two VGT inverse position solutions; given it,r, calculate 
First, however, the gimbal forward solution (Section 4.3.1) is required to start 

the VGT inverse solution. As shown in Fig. 5, the VGT inverse solution is divided into 
solutions A and B, separated by intermediate variables 0 1 ,0 2 ,03. 
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42.1 VGT Inverse A. The cross longeron face angles 0 1 ,02, ©3 are calculated in 
this section given the pointing vector n and the gimbal extension r. Symmetry is exploited 
to obtain the solution. Figure 8 gives the vector diagram for this solution. A similar 
solution was presented by Padmanabhan, et. al., (1992). 



Figure 8 

VGT Inverse A Vector Diagram 


The unknown angles 0i,0 2 ,03 determine the vector locations B Q,,i = 1,2,3, 
respectively (see Fig. 6B. The following three scalar constraint equations are written (see 
Fig. 8): 

a, ■U l = 0,i = 1,2,3 (18) 

Equations 18 state that each B Qj,i = 1,2,3 must lie on a plane parallel to the actuation 

B S 

plane. That is, a, (the vector locating Q in the plane offset by — parallel to the actuation 
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plane) must be perpendicular to U[ (the normal to the actuation plane). From Figs. 7 and 
8, the following vector equation is written: 




where B P — pU j. Equations 18 become: 


B O i -—+-U l \u l = 0,i = 1 , 2 , 3 . 
2 2 


The vectors B Q,,i = 1,2,3 are given symbolically in Eq. 5. Equation 14 is rewritten to 


yield B P : 


B P = r(x B +n). 


r, X B = {1 0 0} T , and n = {n* n y n 2 } T given, B P is known: 


*P = rl n 


Now the actuation plane normal can be calculated: 


"P 

Ul = \\ B Pl 


where: 


B P = rjl+2n x +n x +n y + n 2 =r A /2(l + n^) 


Substituting Eq. 23 into Eq. 20 yields the following three constraint equations. 


*G— £}® ,,+ §M = 0 -'" 1 > 2 - 3 
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Substituting Eqs. 5, 22, and 24, into Eqs. 25, the three constraint equations become: 

EyC\ + + Gj = 0 

E 2 C2 + F 2 S2+G 2 =0 (26) 

^3 C 3 + F3S3 + G3 = 0 


* 


where: 


Ei = Nrn y 
F^Nrll+nJ 


G\ = 


r 

v 


n J ,G lj ,-r(l+n ;( .)+S 



Eo - -Nr\ 


1 

-A.,-- 


V3 


^2 y 2 
F 2 = + 


n T 


G) = r 


n y 0 2y + n z 0 2z - r(l+ n x ) +S 



(26a) 


(26b) 


E-> = -Nr 


1 ^ 

—n v +—n. 


F 3 = Nr{l+n x ) 

n,O h +n,0 3l -r(l+nJ+S^±) 


(26c) 


G3 = r 


Defining «* = 1 + «* , and dividing by r, the three constraint equations simplify to: 


E i c i + Ei s i +G t =0 

E 2 c 2 +F l s 2 +G 2 =0 (27) 

F3C3 + Fis 3 + G 3 = 0 
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(27a) 


E x = Nn y 

Fi = Nn x 

Gi= -^-^ +s 1t 


F 2 =-j(”y-J3”z) 


(27b) 


E 3=~{f'y+J3 n z) 


G ’‘-™’ + T-{jT n y + "‘) +s Jf 


(27c) 


0 1 — r 2t 

Using the tangent half-angle substitution (t = tan—, cos0. = =-, sin0 = -), the 

2 1 +t 2 l+t 2 

general solution to £, cos0, + F i sin 0, +Gj = 0 is given in Eq. 28. 


(0,) u =2tan 


-1 


-F t ±^Ef+lf-Gf 


V , 


(28) 


Due to the multiplicative factor 2 in Eq. 28, the quadrant-specific atan2 function 
need not be used. However, due to the quadratic degree of the tangent half-angle 
substitution, there are two possible solutions for each 0, (the ± in Eq. 28). Since each 0, 

is solved independently of the other two, there are a total of eight (2 3 ) possible solutions 
to the VGT Inverse A problem. This means that in general, eight distinct VGT 
configurations exist to yield the same n,r. For the physical VGT hardware modeled in 
the examples (Section 5), extensive computer simulation reveals there is only one valid 

inverse solution when joint limits are considered. This solution corresponds to choosing 
(0,) 2 ,i = 1,2,3. 

4.2.2 VGT Inverse B. The leg lengths Lj , L2 , L 3 are calculated in this section 
given the cross longeron face angles 0i,0 2 ,03. First, given 0 1 ,0 2 ,03, the vector locations 
B Qi,i =1,2,3 are calculated from Eq. 5. Then the Euclidean norm is used to calculate L,. 
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A=|"a- f 'e 1 l 

^=ir&-"a| 

h = \‘Qi- B Qi\ 


( 29 ) 


Given each 0 1 , 0 2 ,03 from the VGT Inverse A solution, there is a unique 
because only the positive square root is allowed. 


4.3 Extensible Gimbal Position Kinematics 

In order to control the VGT as a gimbal with extension, the kinematic equations of 
this virtual device are developed. Section 4.3.1 presents the gimbal forward solution, 
which is required to start the VGT inverse position solution (see Fig. 5). Section 4.3.2 
presents the gimbal inverse solution, which is required to complete the VGT forward 
position solution. Figure 3 shows the kinematic model for the virtual extensible gimbal. 
The extensible gimbal position solutions relate the input parameters a,J5,r to the output 
parameters n,r (the r transformation is identity). The normal to the top platform is n 
while the gimbal extension is r. The gimbal angles are a,p. The center of the virtual 
gimbal moves with respect to the VGT base. Section 4.3.3 presents the gimbal forward 
transform. Instead of n,r, the output parameters here are [ ^t], the full homogeneous 

transformation matrix (Craig, 1989) relating {P} in terms of {B}. 

4.3.1 Extensible Gimbal Forward. In this section, n,r are calculated given 
a,P,r. Appendix D gives the detailed extensible gimbal kinematic diagram and the 
associated Denavit-Hartenberg parameters. Using these Denavit-Hartenberg parameters 
in the forward kinematics, n is given in Eq. 30. is the unit direction X P expressed 

in {B}: 


h = 



cacp 
jacP > 


-s P 


(30) 


The r transformation is identity. That is, r is viewed here as both an input and output 
parameter for the extensible gimbal. The following equation calculates the value of p, 
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which is the length from the origin of {B} to the origin of {P} (see Figs. 7 or 8). Equation 
31 was derived from Eq. 24. 


p = r > /2(l + cacP) 


(31) 


4.3.2 Gimbal Inverse. In this section, a,P,r are calculated given n,r. Again, 
the r transformation is identity. To solve for a, (5, the given n = {n x n y n z ] is equated 
to the expressions of Eq. 30. A ratio of the y and x terms yields a; cp divides out, but its 
sign is unknown. Therefore, there are two possible a solutions: 


ctj = A tan 2 


/ n \ 

_y_ 

\ n x ) 


(32A 


a 2 = A tan 2 


-n 


\~ n x 


— ' = ^+71 


(32B 


A ratio of the z and x terms yields P ; each a result yields a unique P , so there are two 
possible P solutions: 


Pi 


= A tan 2 


n z ca 

-n r 


(3 3 A 


P 2 =Pi+ti 


(33B 


The primary solution a 1 ,p i is in the range of interest for the VGT. With a 2 ,P 2 , 
the origin of {P} and direction of X P is the same (i.e. a 2 ,P 2 is a valid solution), but Y P 

A, 

and Z P are reversed. 

4.3.3 Gimbal Forward Transform. In this section [pTj is calculated given 

a,p,r. [fr] is the homogeneous transformation matrix (Craig, 1989) giving the position 

and orientation of {P} relative to (B). Since the virtual extensible gimbal has only three 
degrees of freedom, [^r] is a subspace of six-dimensional Cartesian space. The output 

parameters h,r are one possible output set to represent the three VGT degrees of 
freedom, gives the full Cartesian information for the VGT. This transformation 
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matrix is useful for VGT simulation and control. It is also required to calculate the VGT 
node locations (see Appendix E). Given the Denavit-Hartenberg parameters in Appendix 
D, [frj is presented in Eq. 34. Note that the first column contains h = { 5 Xp}. 


c acp 

[*i-5 

0 


-5a 

ca 

0 

0 


carp r(l+cacP) 
5oyP rcacP 

cP -r5P 

0 1 


(34) 


4.4 VGT Node Locations 

Following either VGT Forward or Inverse kinematic position solution, the VGT 
node locations can be calculated. These node locations are required for finite element 
analysis, VGT graphical simulation, and other types of analyses. The node location 
calculations are detailed in Appendix E. 
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5. EXAMPLES 


Three examples are presented in this section to demonstrate the calculation of 
complete forward and inverse position kinematics for the VGT as an extensible gimbal. 
The three examples in this section are for the NASA VGT hardware, shown in Fig. 1 . The 
five nominal parameters for this VGT are (inches): 


Lo = 36 
1 = 34 
l min = 36 
Lmax ~ 55.5 
S = 4.75 

These examples have been verified by computer simulation (inverse kinematics checking 
forward and vice versa) and with hardware measurements. Units for lengths are inches, 
units for angles are degrees, and n and the rotation matrix part of [ fr] are dimensionless. 


5.1 Example 1 

When all input leg lengths are equal, the top platform remains horizontal. All 
variables in Fig. 5 are given below. These numbers may be reproduced with either 
forward or inverse kinematics. 


Lj - Z-2 - Lj = 45 


0j =0 2 = 0 3 = 122.7 

n = {l 0 0} T 
a = P = 0 r = 26.64 


5.2 Example 2 

When two input leg lengths are minimum, and the other is maximum, a maximum 
gimbal angle case results. The possible combinations of L Mm ,L MA % for all three input 
legs lengths yield the extreme gimbal angles (some are in combinations of a, (3). The 
example below gives the maximum positive a for this VGT. 
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Lj =L, = 36 12 = 55.5 


e!=79.1 e 2 = 03 = 138.6 

n = {.673 .739 0} r 
a = 47. 7 p = 0 r = 28.73 


5.3 Example 3 

This is a general example within the motion range of the NASA VGT. For this 
case, [fr] is calculated. Also, Tables 1 and 2 present the multiple solutions of the VGT 

inverse problem. Each set in Table 2 was given as input to the forward 

kinematics, resulting in the same a,(3,r for each solution. As seen in Table 2, only the last 
solution is valid when actuator joint limits 36 ^ ^ 555 are considered. This condition 

has been observed in all of the computer simulations which have been performed, which 
indicates that multiple solutions do not present a problem for control of the VGT as an 
extensible gimbal. 


I, =45 L, =53 La = 50 
0 X =117.9 0 2 =127.4 03 = 141.4 

n = {.955 . 242 .171} r 


a = 14.: 

2 

P = -9.8 

r 

= 24.92 


*.955 

-.246 

-.166 

48.73* 

m- 

.242 

.969 

-.042 

6.04 

.171 

0 

.985 

4.26 


0 

0 

0 

1 


Table 1 Multiple 0, Solutions 



e, 

e 2 

e 3 

1 

47.9 

51.0 

54.2 

2 

117.9 

127.4 

141.4 
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Table 2 Multiple I* Solutions 



A 

h 

h 

1 

14.5 

12.4 

13.5 

2 

14.5 

30.2 

29.7 

3 

24.7 

25.3 

13.5 

4 

24.7 

53.0 

29.7 

5 

21.3 

12.4 

21.5 

6 

21.3 

30.2 

50.0 

7 

45.0 

25.3 

21.5 

8 

45.0 

53.0 

50.0 


6. CONCLUSION 

This paper presents the complete forward and inverse kinematic solutions for 
control of a three DOF double octahedral VGT module as an extensible gimbal. The 
results can be used for design, analysis, simulation, and control. A VGT can be used as a 
joint in a large, lightweight, high load-bearing VGT-based manipulator for earth- and 
space-based remote operations, plus industrial applications. 

Using the Stewart platform analogy, there are eight possible solutions to the VGT 
forward position problem. However, in this paper, a numerical method is used for 
solution, which yields one solution, closest to the previous configuration. There are eight 
solutions to the VGT inverse position problem. All eight solutions are calculated in this 
paper, but it was observed that only one solution is valid when practical joint limits are 
considered. 

. The results of this work have been used to control the NASA VGT hardware as an 
extensible gimbal, to demonstrate the capability for this device as a joint in a VGT-based 
manipulator. This work is an integral part of a comprehensive VGT-based manipulator 
design, simulation, and control tool. 
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APPENDIX A. VC.T Mobility 


The Kutzbach equation (Mabie and Reinholtz, 1987) is useful for calculating the 
spatial mobility (Af) for a given device: 


M = 6{N-l)~ 5 J t - 4 J 2 - 3 / 3 - 2/4 - J 5 (A. 1 ) 

where M is the mobility , or number of degrees of freedom of a device, N is the number of 
links including the fixed link, and /,• is the number of /-DOF joints. 

For a theoretical double octahedral VGT, all joints are modeled as spherical joints, 
and any strut may be linearly actuated. For the case of three actuators on the mid-plane, 
the mobility is: 


M = 6(20-1) -5(3) -3(27) =18. (A.2) 

There are twenty links (top and bottom plates, twelve cross longerons, and two for each 
actuated link), three prismatic joints, and twenty-seven spherical joints (for each mid-plane 
node, six links theoretically join at one point; so five spherical joints are required). 
However, each link with spherical joints on each end has an idle DOF, rotation about its 
own axis. There are fifteen idle freedoms (twelve cross longerons and three actuators), so 
the effective mobility is M = 18-15 = 3. 

For the specific NASA VGT hardware, the design goal was to replace as many 
spherical joints as possible while maintaining three DOF. To enumerate links and joints, 
consider a VGT subchain: one pair of cross longeron faces connecting the base plate to 
the top plate; e.g. connecting B 3 B X to B 3R B lR through Q 3 and Q 3R (see Fig. A.l). As 

shown in Fig. 6 b, two cross longerons forming a face may be replaced with a kinematically 
equivalent single link, rotating about the base (or top) plate with a single revolute joint. 
The two cross longeron faces are connected by two universal joints, separated by the 
yoke. An idler rotates relative to the yoke. The two actuated links connect to the idler via 
two revolute joints. Therefore for each such subchain, there are six links, two universal 
joints (which are two-DOF joints), and five revolute joints. For the entire VGT, there are 
twenty links (six links per subchain times three subchains, plus the top and bottom plates), 
six universal joints, and eighteen one-DOF joints (five revolute joints per subchain times 
three subchains, plus three actuated prismatic joints). Therefore: 

M = 6(20 - 1)-5(18)- 4(6) = 0 (A.3) 
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This result predicts that the NASA VGT is a statically determinate truss with no freedom 
to move! The Kutzbach equation is known to fail when the device contains special planar 
geometry within its spatial arrangement. The VGT actuation mid-plane is a triangle. The 
actuator links connect to the idlers, which are always perpendicular to the actuation 
triangle. The Kutzbach equation does not know of this special arrangement and correctly 
predicts zero DOF (if these three idler axes were skew to each other, no motion would be 
possible). To get the correct results in the Kutzbach equation, the one-DOF prismatic 
actuators must be modeled as two-DOF cylindrical joints. This removes three joints from 
J x and adds them to J 2 : 


M = 6(20-1) -5(1 5) -4(9) = 3 (A.4) 

Thus the NASA VGT has three DOF. 



Figure A.l VGT Subchain 
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APPENDIX B. 
Fixed VGT Geometry 


This appendix presents fixed VGT geometry for use in forward and inverse 
position analyses. Referring to Fig. 6b, the points B O t = {o O iy O i2 } are the virtual 

revolute joint locations in the Stewart platform analogy. Fig. B.l shows the details 
necessary to calculate these locations. 



0 


0 


0 

-v 

II 

Lo 

2V3 

0 

b 0 2 = < 

Lo 

4>/3 

Lo 

. 4 j 

Bn _ 
0 3 =■ 

V 

4V3 

io 

4 



<|)i =180 
=300 
= 60 


As shown in Fig. B.2, the top platform geometry is identical to the base. 
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Figure B.3 shows one of the cross longeron planes. These three planes are 
identical, so: 

M-w-fiZ. (b.2) 

Qi 



Lower Cross Longeron Face Geometry 
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Figure B.4 shows that the upper cross longeron planes are symmetric to the lower. 



Qir 

Figure B.4 

Upper Cross Longeron Face Geometry 
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This appendix presents numerical solution of Eqs. 7 using the Newton-Raphson 
method. These equations, arising in the VGT Forward A position analysis, are coupled 
and non-linear (transcendental). The general Newton-Raphson method for numerically 
solving non-linear systems of equations is discussed in Dahlquist and Bjorck (1974, 
Section 6.9.2). The Newton-Raphson method is applied by Williams (1992, Appendix C) 
to solve the forward position problem of the Stewart platform. In this appendix, the 
Newton-Raphson method is adapted to solve the VGT Forward A problem. It requires an 
initial guess which is "sufficiently close"; quadratic convergence is guaranteed under this 
condition. In real-time control, an excellent initial guess is the previous known 
configuration of the VGT. The VGT has eight possible forward kinematic solutions, but 
the method in this appendix only tracks one, closest to the previous configuration. The 
inverse Jacobian matrix is presented symbolically, so Gaussian elimination is not required. 

Equations 7 are three coupled equations in transcendentals of the unknown cross 
longeron face angles © = {0! 0 2 0 3 } 7 ', rewritten in Eqs. C.l and C.2. 

Jv(0)=O;/ = 1,2,3 (C.l) 

F x (0) = q + c 2 + AciC 2 - 2 Asfa + B = 0 

F 2 (0) = c 2 + c 3 + Ac 2 c 3 - 2 As 2 s 3 + C = 0 (C.2) 

F 3 (0) = c 3 + q + Ac 3 q - 2As-i,Si + D = 0 

The general Newton-Raphson iteration for step k+1 is summarized below. 

JdQ k =-Fi{e k ) (C.3) 


The Jacobian matrix is a multi-dimensional form of the derivative: 



(C.4) 


Equation C.3 is solved for 5©*, then the updated vector of unknowns is calculated with 
Eq. C.5. 
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(C.5) 


0 *+i =0 * +50 Jt 


Iteration continues until |5@*fl<6, a user-defined tolerance. The method requires an 
initial guess, @ 0 , "sufficiently close" for convergence. 

For the VGT, the Newton-Raphson Jacobian matrix is given below. 


*1 

* 3 

0 

0 

k 4 

*6 

*2 

0 

*5 


(C. 6 ) 


where: 

*1 = KqSi + K^Ci 

K 2 — K c s x + KpCi 

*3 = %A S 2 + *Z) C 2 
*4 = * 0*2 + * F C 2 
*5 = *A 5 3 + *D C 3 
*6 = *B^3 + KgCs 


K a = -{1 + A Ci ) 
tf B = -(l + Ac 2 ) 
/sT c = -(l + Ac 3 ) 

Kp — ~~2>jA.Si 

Kg — — 2As 2 
Kg — —2As$ 


This Jacobian matrix is inverted symbolically, so Eq.C.3 is solved analytically as follows: 

60*=-y _1 F;(0*) (C.7) 


where: 



DET 


*4*5 

— *3*5 

*3*6 

*2*6 

*1*5 

-*1*6 

— *2*4 

*2*3 

*1*4 


(C. 8 ) 


det=k 1 k 4 k 5 +k 2 k^k 6 


(C.9) 
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APPENDIX D. Virtual Extensible Gimbal Kinematic Details 


This appendix presents the Denavit-Hartenberg parameters and the associated 
kinematic figure for the virtual extensible gimbal. This information was used to derive 
for Sections 4.3.1 and 4.3.2, and [pT] in Section 4.3.3 and Appendix E. For the 

coordinate frame assignments shown in Fig. D.l, Table D.l gives the Denavit-Hartenberg 
parameters (Craig, 1989, convention) for the extensible gimbal. The gimbal angles are 
a,P and r is die symmetric gimbal extension. 

Table D.l Extensible Gimbal D-H Parameters 



Figure D.l Extensible Gimbal Kinematic Details 
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APPENDIX E. 
VGT Node Locations 


Following either the VGT forward or inverse kinematic position solution, the VGT 
node locations can be calculated. These node locations are required for finite element 
analysis and VGT graphical simulation. Figure E.l shows the 12 VGT nodes calculated in 
this appendix. There are four sets of nodes: base nodes fi, , lower actuation plane nodes 
, upper actuation plane nodes Q R , and top platform nodes B iR (t = 1,2,3) respectively. 

a 

this appendix these nodes are calculated given 0 1 , 0 2 , 03 , I/j, ot,p,r. All nodes are 
expressed with respect to the fixed {B} frame. 
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E.l Base nodes £,;/ = 1,2,3 

These are fixed for a single VGT cell. The geometry is shown in Fig. B.l. 


r 



0 


o i 

< 2V3 

II 

~Lo 

iS 

A) 


-A) 

, 2 


2 J 



(E.1) 


E.2 Lower actuation plane nodes Q { ;i = 1,2,3 

These expressions were required for the VGT Forward A and VGT Inverse A 
solutions in Sections 4.1.1 and 4.2.1, respectively. Equation 5 is repeated in Eq. E.2. 
Appendix B presents fixed VGT geometric values used to calculate these nodes. 
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E.3 Upper actuation plane nodes Q iR ,i = 1,2,3 

These three nodes are Q^i = 1,2,3 reflected about the actuation plane of 
symmetry. They are calculated by extending the VGT joint offset along the actuation 

A 

plane normal U l . 

B Q,R= B Q, +SU t (E.3) 


E.4 Top platform nodes B iR’i -1,2,3 

These three nodes are B,;i = 1,2,3 reflected about the actuation plane of 
symmetry. They are calculated by transforming the top platform nodes from the {P} frame 
to the {B} frame. 


Bn BrpPn 

BiR-pl iS iR 


(E.4) 
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The transformation matrix %T, developed in Section 4.3.3, is repeated below. 



cacp 

-s a 

carp 

r(l+cacP) 

sacp 

c a 

jcwP 

rjacP 


0 

c P 

-rs p 

0 

0 

0 

1 


(E.5) 


From symmetry, the top platform nodes in the {P} frame are identical to the base nodes in 
the {B} frame (see Figs. B.l and B.2). 

P BiR= B Bi (E.6) 
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